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Abstract. The mcsoscopic concept is applied to the description of microc- 
racked brittle materials. The mesoscopic equations are solved in a special case 
when the microcracks are developing according to the Rice-Griffith evolution 
law. The evolution of the crack distribution function is investigated in case 
of simple loading conditions and for two different initial crack distribution 
functions. The time dependence of the average crack length is calculated, too. 



1. Introduction 

There are several different approaches for the description of the behavior of ma- 
terial with complex microstructure, e.g. granular media, liquid crystals or damaged 
solid materials. One of the basic approaches is the macroscopic one applying ther- 
modynamic concepts and principles to describe the microstructure (e.g. internal 
thermodynamic variables). The other basic but different approach is the statistical 
one where the microstructure is considered as a statistical ensemble and the macro- 
scopic variables are calculated from the interactions between the microstructural 
elements using the principles of statistical mechanics. Modelling the microstructure 
tests the limits of applicability of both theories and therefore the basic principles 
of our understanding regarding the material behavior. 

The microscopic elements of statistical cracking phenomena are the single cracks. 
The form and validity of single crack evolution equations is an important problem 
in statistical microcracking. Beyond the well known classical results (see e.g. 
there are new and promising approaches and investigations in this field that should 
be considered. Langer and his coworkers EJ E] gave a microscopic background 
of the internal variable approach of plasticity ^k. Using molecular dynamic sim- 
ulations they found some good microscopic candidates of the internal variables. 
Based on these observations a so called shear deformation zone (STZ) model was 
established, suggesting a dynamical equation of these internal variables, including 
also special unilateral requirements. Several inconsistencies of classical single crack 
models can be explained using this theory. Another interesting approach to single 
crack phenomena considers cracking as a kind of nucleation problem, where elastic 
and conservation effects are included in the nonlocal (coarse grained) phase field 
model El Hj , or the elastic effects are included into a more direct nucleation 
approach ^JJ ^2 1X3 • Both STZ and the phase field approaches can be considered 
as phenomenological thermodynamic theories. 

On the other hand, partially independently of the previous considerations micro- 
cracking is treated also from a statistical point of view. There are two competing 
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models in this respect. Sornette and his coworkers considers cracking as a crit- 
ical phenomenon where rupture is a critical point. Their approach is supported 
by experimental scaling observations, by numerical simulations (see [131 114| and 
the references therein). Other observations and simulations suggest that the rup- 
ture of the material is better considered as a first order phase transition and the 
cracking as spinodal nucleation |15| . Scaling calculations based on the previously 
mentioned nucleation models partially support that approach. From the point of 
view of nonequilibrium thermodynamics the background of both the statistical and 
single crack models is the violation of thermodynamic stability in the nonequilib- 
rium state space. This observation alone can unify several classical failure and 
fracture criteria and explain experimental observations 117) . 

However, between pure phenomenological and statistical approaches there is an 
independent family of theoretical descriptions that we can call hypercontinuum the- 
ories. In these theories the physical quantities are not simple macroscopic fields 
given on the appropriate space-time (classical, special or general relativistic), but 
one assumes that in addition to the space-time variables they depend also on some 
mesoscopic variables characterizing the microstructure of materials. These theo- 
ries have almost as long history as the traditional macroscopic and microscopic 
ones. At the beginning they were used to model rotating internal particles in mi- 
crocrystalline materials (see the basic works of the Cosserat brothers and Mindlin 
|181 119j). nowadays generalized continua are used to model e.g. granular media 
[231 EH 122, liquid crystals [212 El, ferrofluids [2H1EE], polymer solutions [23EHj, 
but one can meet hypercontinuum in general relativity, too [2S EOj (see also [S3 
and the references therein) . 

Hypercontinuum theories represent models of the material where one can get 
information from below the macroscopic but above the microscopic level. That gives 
the basic motivation and modelling idea to the so called mesoscopic theory, that uses 
both macroscopic, continuum and microscopic, statistical concepts constructing the 
basic equations of motion for the hypercontinuum. This approach was successfully 
applied for the description of liquid crystals (see [32] an d the references therein) . 

Recently a mesoscopic theory was suggested for a special class of damaged mate- 
rials, for microcracked brittle solid media |33| . In this paper we will see, that consid- 
ering the peculiar properties of microcracking simplifies the mesoscopic equations 
well enough to make them solvable with relatively little effort. This simplification 
results in a generalized Liouville equation that we will solve considering a simple 
evolution equation for the microcracks in case of two different initial crack distri- 
bution functions. We applied Rice-Griffith microscopic dynamic equations as the 
simplest unilateral microdynamics that is exactly (analytically) solvable on the sta- 
tistical level. The explicit expression of the microcrack distribution function makes 
possible to investigate the relevance of some damage parameter and its respective 
evolution equation and the role of microscopic unilateral constraints. 

2. Basic concepts of mesoscopic microcracking 

In this section we give a short summary of the mesoscopic microcrack theory and 
the related equations governing the development of crack distribution function. For 
a more detailed treatment we refer to |33| . 

In a mesoscopic theory of a microcracked continuum our quantities are inter- 
preted on the directional space that includes the 1 length of the microcracks and 
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the x, t space-time coordinates. The cracks are considered as penny shaped surfaces 
with area irl 2 and orientation n. These two quantities characterizing the crack are 
put together as a vector 1 = In, where I is the crack radius. The microstructure con- 
sists of this surface vector of the cracks. This vector is an axial one, that we do not 
consider as an antisymmetric tensor but for the sake of convenience we use normal 
polar vectors with a symmetry requirement on their number density function: 

W(l,x,t) = W(-l,x,f). 

The basic equations of the mesoscopic theory are the mesoscopic balances of the 
fundamental physical quantities. Due to the enlarged configurational space, the 
mesoscopic velocities and the mesoscopic material current densities have six com- 
ponents and can be decomposed into a 'normal' part (with respect to the position 
variable) and a 'mesoscopic' or 'directional' part. The macroscopic number density 
of microcracks of any size and orientation at position x and time t. N can be in- 
troduced with an integration over the directional part of the mesoscopic space |H4| : 
N(x,t) := (JV(l,x,t)) := 5 L 3 iV(l, x, t)dVt. Now the balance of the mesoscopic 
number density, the generalized continuity equation can be written as: 

(1) ^ + V -(iVv) + Vr(iVvO = a„. 

Here v is the translational velocity, the space part of the mesoscopic velocity and 
v; is the 'directional' part. V; denotes a derivation with respect to the variable 1 
and as a part of the mesoscopic divergence represents physical requirement that the 
number of crack with a given length at a given material element can be changed by 
changing the length of some smaller or larger cracks. • dot stands for the divergence, 
the source term a n characterizes the creation and coalescence of microcracks. 

We introduce a distribution function / as the probability density of a crack 
having a particular length and orientation as follows 

2 f(l, n, x, t) = - 

N(x.,t) 

From the mesoscopic balance of number density follows the equation determining 
the time development of the crack distribution function is 

(3) ^ + V • (/v) + V ; • (/v«) = ^.(l +v .vjN + a n 

In addition we have balance equations for the mesoscopic densities of all other 
extensive quantities e.g. one can give the mesoscopic balances of mass, momentum, 
angular momentum and energy |85| . The macroscopic balances can be introduced 
with averaging the corresponding physical quantities over the additional, enlarged 
part of their domain, using the previously introduced directional probability distri- 
bution function. 

We are interested here in the evolution of crack lengths and damage is described 
in terms of the crack length distribution function. However, crack orientation was 
introduced as an additional mesoscopic variable together with crack length, and 
there may be an arbitrary distribution of crack orientations. Special cases of ori- 
entation distributions are: an isotropic distribution of all orientations, or, as the 
other extreme case parallel orientations of all cracks. Very often we will meet cases 
in between: the cracks are not oriented exactly parallel, but there is a preferred 
orientation. In these cases the orientation distribution can be rotation symmetric 
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(uniaxial) or not (biaxial). This classification in terms of the orientation distribu- 
tion is analogous to liquid crystals |36l I37| , where the mesoscopic variable is the 
orientation of elongated molecules. Here we want to consider only averages over 
all crack orientations. We introduce the crack length distribution function as the 
average over all crack orientations: 

(4) f(l,x,t) = ±- [ f(l,x,t)d 2 n . 

4?r J S 2 

Averaging of the differential equation ||3J) over orientations results in the equation 
of motion for the crack length distribution (|llfl . (For more details see It is 

the same form of differential equation for any orientation distribution. Only the 
value of the crack length change velocity vi depends on the orientation distribution 
via the effective stress acting on cracks of different orientations. Therefore we have 
only a parametric dependence of the differential equation to be solved for the crack 
length distribution on the orientation distribution, and the same form of solutions 
is expected for any orientation distribution. 
The distribution functions are normalized: 

(5) / / f{l,x,t)l 2 dld 2 n = l Vx,t 
Js 2 Jo 

and 

/>oo 

(6) / f(l,x,t)l 2 dl = l Vx,t . 
Jo 

3. Solution of the mesoscopic equations 

There are essentially two methods to gain information from the basic equations 
of the mesoscopic theory. The first one, that was used mainly for liquid crystals, 
is to calculate a moment series expansion of the equations. Using this expansion 
we can introduce a hierarchy of macroscopic internal variables together with their 
evolution equations. This can be useful also in our special case, the nature of 
macroscopic damage parameters and the properties of their evolution equations is 
an old and strongly discussed problem in continuum damage mechanics. However, 
now we will prefer another method: we will solve the mesoscopic equations in simple 
cases to get direct information how the distribution function of microcracks evolves 
in time. 

The system of mesoscopic balance equations given in the previous section is too 
general for most of the practical problems and materials in damage mechanics. In 
case of brittle (rock or ceramic) material the following simplifications are suitable 

(1) The base material does not have an internal spin, that is a crack does not 
rotate independently from the base material. 

(2) There are no couple forces and couple stresses. 

(3) There are no external body forces (f = 0), 

(4) The material is in (local) mechanical equilibrium (v = 0). 

(5) The translational velocity v does not depend on the crack size and orienta- 
tion. It is the same for all cracks, namely it is equal to the velocity of the 
center of mass of the surrounding continuum element (v(l, x, t) = v(x,t)). 
This does not imply restrictions on the length change velocity. 

(6) There is no creation of cracks. This is not a crucial restriction. For example 
preexisting voids are considered as microcracks of very small size. In the 
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progress of damage they are growing, but their number is constant. (The 
voids in the initial stage are "counted" as cracks.) 
(7) There is no crack coalescence. Together with the previous assumption we 
postulate that there are no crack sources, i.e. the number of cracks is 
constant in time. Let us remark, that the absence of coalescence does 
not mean the absence of interaction. Interaction can be considered in the 
microscopic dynamics as we will see below. 
Let us remark that these assumptions are not necessary, more general cases can 
be treated in the frame of our model. E.g. crack creation and coalescence can 
be considered through a suitable source term in l|T]l. However, in this paper we 
maintain the theory as simple as possible to investigate some effects of microscopic 
unilateral dynamics. 

All the previous assumptions are reasonable in case of solid materials. Consider- 
ing these assumptions we find that the balance of moment of momentum is trivial, 
the balance of momentum simplifies to the condition of mechanical equilibrium: 

V-t = 0, 

where t is the stress. Moreover at present we are not interested in the changes of 
the mass density and the produced energy during cracking, we will concentrate on 
the evolution of the crack distribution function. 

As a result of the conditions 5-7 the macroscopic number density of microcracks 
does not change in a material element (only the lengths can be changed) and there 
is no source term at the mesoscopic level. Therefore there is no source term in 
the differential equation for the crack distribution function J3J). Due to condition 4 
choosing a comoving frame results in a zero translational velocity and we get the 
following simple equation: 

(7) §£ + Vr(/v*) = 0. 

Let us mention here that in this case the equation above has a direct interpre- 
tation. From a pure statistical mechanical point of view this is a continuity equa- 
tion for the probability density of microcracks having the " microscopic" dynamical 
equation 

i = v,(i). 

Here the point denotes the substantial (comoving) time derivative. However, the 
special form of the differential equation for the distribution function depends on the 
equation for the length change velocity. This expression can involve a dependence 
on orientations and lengths of surrounding cracks and also several other forms of 
crack interaction. The equation JJJ is a generalized form of the classical reversible 
Liouville equation of statistical mechanics, determining the dynamics of the distri- 
bution function / with the help of a known deterministic (!) microscopic dynamics 
|39| . From this observation it is clearly seen, that in the simplest non-interacting, 
reversible case the mesoscopic theory is fully compatible with the basic statistical 
theories and continuum theories, too. 

4. Rice-Griffith dynamics 

As an example we will use one of the simplest, but well known dynamics that 
was given by Rice from general thermodynamic considerations giving an irreversible 
thermodynamic background to the energetic considerations of Griffith The 



6 



P. VAN 1 , C. PAPENFUSS 2 , W. MUSCHIK 3 



condition is derived for a single crack under uniform tensile loading parallel to the 
direction of the crack surface normal. The main assumption in the construction of 
the dynamics is that the Gibbs potential depends not only on the stress but also 
on the crack length in the following way 

(8) G(a,l) = al-^a 2 l 2 -la 2 . 

Here a is usually considered as a kind of surface energy, j3 is connected to the 
energy release rate in case of simple geometric situations, 7 is the elastic coefficient, 
I is the length of the crack and a is the tensile stress parallel to the crack normal. 
Let us remark here that similar conditions can be given in case of compressive 
stresses with modified material parameters (a compressive stress can result in a 
tensile force on a crack surface for certain crack orientations) . The first term in © 
is usually interpreted a surface energy of the crack, the second as the extra elastic 
energy of the surrounding continuum and the last one is the pure elastic term in 
this simple case. The sign and Legendre transformation conventions are according 
to Landau and Lifsic |42j . 

The parameters a, (3 and 7, the coefficients in the Gibbs potential, are effective 
parameters of the surrounding continuum, including other cracks. 

From this Gibbs potential we can derive the entropy production a s due to the 
crack extension in case of isothermal conditions 

d 2 G\~ 1 fdG \ . dG- 



da 2 J \da J dl 

It is easy to identify thermodynamic forces and currents in the expression. More- 
over, we are interested only in the crack propagation effects, therefore we will sup- 
pose that during the propagation of the crack the material is in mechanical equi- 
librium, because the crack propagation is considerably slower than the deformation 
process. We conclude, that 

(9) e=- =(7 + 

OCT 

where e is the deformation and we arrive to the evolution law of crack propagation 

l = -L^- = -a' + P'la 2 . 
al 

Here L is a positive kinetic coefficient (according to the Second Law) and a' = 
La, P' — Lj3. The equilibrium solution of this equation gives the celebrated Griffith 
condition 40 . After simple stability considerations one can see easily that the 
cracks start to grow if the right hand side of the above differential equation is 
positive (at least if we are investigating brittle rocks under room temperature) 
expressing the unilateral nature of crack propagation. If the healing of the cracks 
is excluded, then the above evolution law can be written more exactly as 



(10) I 




P 'la 2 if a' < (3' la 2 , 
otherwise. 



This evolution equation will be a core for the mesoscopic dynamical investigation 
resulting in an explicit expression for the mesoscopic velocity. The above unilateral 
dynamics will be called Rice-Griffith dynamics. 
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Let us remark here that there is much more information that can be gained 
from these thermodynamic considerations. For example as a condition of ther- 
modynamic stability we get an upper limit for the crack length > £ 2 ) that 
one cannot read from the single Hl()|) evolution equation without considering the 
mechanical interaction. The role of thermodynamic stability in the failure of real 
materials is investigated in more detail in [1611171 l4"5] . where it is shown that local- 
ization and critical damage can be considered as the loss thermodynamic stability, 
and some empirical failure criteria (like Griffith) can be understood well from a 
thermodynamic point of view. 



5. Solutions of the Liouville-approximation 

Let us consider now a sample with tensile hydrostatic loading. The mesoscopic 
variable in a given frame can be decomposed into the crack length and the crack 
orientation 1 = (7,0). Similarly the mesoscopic velocity can be decomposed into 
a length change part and an orientation change part, as v; = (vi,uj). According 
to assumption 1 and 4 the cracks cannot move independently of the surround- 
ing continuum, and let be the material at rest v = 0, therefore the orientational 
change velocity is zero (u> = 0). Here we will investigate pure tensile loading con- 
ditions, when the crack surface normal are parallel to the applied loading. Due to 
our simplifications a the role of crack orientation changes can be neglected in the 
development of the different order parameters [HHj • 

In the following we will investigate constant loading rate conditions a — v a t 1 
v a = const, common in standard failure tests. Considering all of our assumptions 
regarding the cracks we can conclude that the crack density distribution depends 
only on time and the length of the cracks. The evolution equation Q in spherical 
coordinates and in the mesoscopic space is the following 

ni s d fun- 1 dl 2 vi{l,t)f{t,l) 

( ' dt [,) ~ P dl 

where the mesoscopic speed, the crack length change velocity vi is given by the help 
of the right hand side of the equation i|l(J|) . that is U; = I. 

Easy to check that the general solution of the first order partial differential 
equation l|ll|) with the microscopic dynamics i|10|l is the following 

^^^[^-^^(ie-^VT^^^rClAO,^ 3 )) if a'<P'vllt\ 

[/(7,0) otherwise. 

Here T is the generalized incomplete Gamma function defined by T(a, 6, c) = 
Fj(a, c) — Ti(b, c), where Yi denotes the incomplete Gamma function. F is an ar- 
bitrary function to be determined from the initial conditions. One can observe 
that there is no stationary solution of l|ll|) . but the crack distribution with a fi- 
nite support does not change under a definite loading. A solution with arbitrary 
initial conditions cannot be given analytically because of the non-healing micro- 
cracks, the non differentiable right hand side of (|10fl . As a consequence of this 
non-differentiability, the solutions of the partial differential equation can be split- 
ted into two parts connected by a moving boundary. Smaller and smaller cracks 
start to grow arriving at the moving boundary. The arbitrary function F should be 
determined in such a way that the final distribution function will be continuous. 
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Figure 1. The crack length distribution function at different in- 
stants. Exponential initial condition. 

In the following we introduce dimensionless variables with measuring the time 
in units of i = zv 2 ^ 3 , the length in 1 = av a 2 / 3 p~ 1 / 3 anc l the mass m = 
av a °/ 3 /3 -4 / 3 as the crack length I, the stress a, the stress change velocity v a and 
the material parameters a 1 and (3' are measured in SI units. The dimensionless 
microdynamic equation l|10(l does not contain material parameters (e.g. can be 
written as v a = 1, a' = 1 and j3' = 1 in (JTOJ). In the following all the variables are 
understood as dimensionless. 

The solution of this equation is calculated numerically in case of two different 
initial length distribution functions. First an exponential initial crack length dis- 
tribution function was considered 

You can see the solution with the parameter value S = 1 at different instants on 
figure n O n the vertical axis the function l 2 f(l,t) is drawn (because / multiplied 
with the spherical weight has the direct probability interpretation). One can ob- 
serve, that the breaking point of the initial exponential distribution moves to the 
left, according to Rice-Griffith dynamics. 

In figure [21 we can observe the full time and crack length development of the 
distribution function. The projection of the thick pointed line, where the surface 
is broken, onto the I — t plane gives the Griffith condition, considering the time 
dependence of stress. The other lines show the distribution function at the instants 
on figure 2. 
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Figure 2. The time development of the crack length probability 
distribution. Exponential initial condition. 



As a more realistic example, that considers an upper limit of the possible crack 
lengths we treat a bounded uniform initial crack distribution function, too: 

I otherwise. 

In the following we choose l c = 1. In figure [3] one can see the crack length 
distribution at different instants. Figure 0] shows together the time and length 
dependence. The pointed line and its projection onto the I — t plane is the Griffith 
condition, that is more apparent on the backward view. The thin lines are the crack 
distribution functions at the same instants as in figure 4. 

Here the resulted solution is discontinuous, due to the discontinuous initial con- 
dition. The initial discontinuity, the largest cracks, move to larger crack sizes with 
increasing time. 

In |38j we suggest several possible macroscopic internal variables from a meso- 
scopic point of view. Different order and damage parameters are introduced and 
the possible advantages and disadvantages are discussed. Here we will investigate 
one of the most remarkable candidates, the average crack length and calculate its 
time dependence in case of the treated two initial crack length distributions. The 
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Figure 3. The crack length distribution function at different in- 
stants. Stepwise initial condition. 





Figure 4. Different views of the time development of the crack 
length probability distribution. Stepwise initial condition. 



average crack length will be called in the following as damage and defined as 

D(t)=<l>= lf(l,t)l 2 dl= l 3 f(l,t)dl. 
Jo Jo 

The notion of an order parameter is chosen in analogy to the mesoscopic theory 
of liquid crystals, where the analogously defined orientational moments of the dis- 
tribution function, the alignment tensors, are the order parameters in a Ginzburg- 
Landau theory. There the gradient of the order parameter contributes to the free 
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energy density in the form of the Frank elastic energy. Here, in the mesoscopic 
theory of damage, we did not take into account any spatial gradients at all, and 
therefore there is no direct comparison to a Ginzburg-Landau theory. Definitely 
the damage parameter introduced here is different from the order parameter, the 
displacement field in phase field models. The spatial structure of the crack is out 
of the scope of the mesoscopic theory, and we are dealing only with microcracks 
smaller than the typical dimension of the continuum element. On the other hand 
the onset of damage has been interpreted as a loss of thermodynamic stability with 
the help of a free energy depending on a damage parameter density |16| . The 
inclusion of spatial inhomogeneities of the damage parameter, and the formulation 
of an extension of a Ginzburg-Landau type theory of damage (with a mesoscopic 
background) is an outstanding problem for the future. 

After calculating the time dependence of the damage (the dots on figure and H3) 
we can observe that time development is similar to the solution of the microscopic 
dynamic equation. Therefore we can investigate, whether the damage development 
is determined by a Rice-Griffith dynamical equation. More properly we suppose 
that the evolution equation that determines the damage dynamics is the following 

(12) D(t) = -a + bt 2 D{t). 

One can get these equations by averaging the mesoscopic dynamic equations 
without considering the unilateral property of the microdynamics. The above equa- 
tion is the only reasonable candidate to describe the macroscopic dynamics of the 
averaged crack length. 

To compare the exact and the suggested macrodynamics we fitted the solutions 
of the above equation with the exact solution to determine the parameters of the 
damage dynamics and the dependence on the initial conditions. In both cases we 
considered the possibility, that the damage starts from the initial value determined 
from the initial crack distribution (dashed lines on figure and or the initial 
damage belongs to the parameters to be determined (continuous lines on figure 
and|BJ). The results are summarized in the table. 





to 


D(t ) 


a 


b 


Exponential 





1 


0.395 


1 







1.2226 


0.5681 


1 


Stepwise 


1 


0.5 


0.7006 


1 




1 


0.5534 


0.8061 


1 



The parameters were determined by an asymptotic approximation method con- 
sidering more and more sample points. We can observe that the /3' = 1 parameter 
of the microscopic crack dynamics is the same as the analogous b in case of damage 
dynamics, however the a = 1 and a parameters (the "surface energy" times the 
dynamic coefficient) are smaller for both initial distributions. The averaged macro- 
scopic equations especially in case of fitted initial conditions give remarkable well 
approximation for larger times as one can see e.g. on[S] However, at the beginning 
when the number of "frozen" cracks is large, the damage developments are quali- 
tatively different (see figure 0. One can observe a kind of "pseudo healing", the 
damage curves calculated from the macroscopic considerations start to decrease, 
instead of the expected increasing. This effect is due to the unilateral micrody- 
namics. The averaged distribution function is monotonously increasing but the 
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0.5 1 1.5 2 2.5 3 

Figure 5. Time dependence of damage, in case of exponential and 
stepwise initial distribution. The dots are related to the averaged 
distribution function, the dashed and continuous lines are the lines 
calculated from the macrodynamic equations. 



Da 




Figure 6. Time dependence of damage, in case of exponential 
and stepwise initial distribution. (Enlarged initial part of figure[3J) 



suggested approximate macrodynamic equation at the beginning start to decrease. 
The decreasing depends on the initial distribution. The initial part of the curves 
related to the stepwise distribution is an extrapolation in both figures to show that 
an unilateral macrodynamic equation would not improve the situation. 
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6. Discussion and final remarks 

In this paper the equations of mesoscopic microcracking were solved in case of 
simple cracking processes. We have seen that under reasonable assumptions the 
seemingly complex mesoscopic balances are reduced to a solvable first order partial 
differential equation giving a kind of Liouville-approximation of the mesoscopic 
theory. There we could get the mesoscopic velocities introducing some constitutive 
assumptions based on known single microcrack dynamics. The partial differential 
equation has an independent statistical physical interpretation and can be useful to 
describe the evolution of microstructure without the underlying mesoscopic theory, 
too. Moreover, it gives a sure basis of the internal variable calculations, giving 
particular interpretations of the different internal variables and clues regarding 
their dynamics, the changes of the elastic moduli or the stress distribution without 
any further physical assumptions. 

The development of the average crack length was investigated in detail and we 
have got a macroscopic dynamics that was almost the same as for single cracks 
with the remarkable property that the corresponding parameters of the damage 
evolution equation strongly depend on the initial crack distribution. Moreover, 
the statistical calculations resulted in a new qualitative phenomena: the average 
crack length was decreasing at the beginning, seemingly contradicting to the uni- 
lateral microdynamics of the microcracks. That calculations give some insight into 
the applicability of internal variable approaches in noncquilibrium situations. The 
unilateral microscopic conditions alone (without quenched disorder) resulted in a 
macroscopic dynamics of the damage (internal variables) that proved to be a good 
approximation of to the microscopic dynamics, but with parameters strongly de- 
pending on the initial crack distribution. 

Our assumptions regarding the microscopic, single crack propagation evolution 
equations can be extended to far more difficult cases, for example we can con- 
sider more general loading conditions, a microscopic dynamics with terminal crack 
velocity, too. 
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